Dirac monopoles and dipoles in ferromagnetic spinor Bose-Einstein condensates 
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We investigate a radial spin hedgehog, analogous to the Dirac monopole, in an optically trapped 
atomic spin-1 Bose-Einstein condensate. By joining together a monopole-antimonopole pair, we 
may form a vortex line with free ends. We numerically simulate the three-dimensional dynamics 
and imaginary time relaxation of these structures to nonsingular textures and show they can be 
observable for realistic experimental parameters. 
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INTRODUCTION 

Topologically interesting structures in quantum fields 
range from vortices and vortex lattices in single compo- 
nent fields , to Skyrmions in multi-component fields • 
They have long been important in the study of superfluid 
physics and quantum cosmology 0. Experimental 
dilute gas Bose-Einstein condensates (BECs) allow these 
structures to be investigated with unprecedented flexibil- 
ity. 

The successful trapping of atomic BECs by purely op- 
tical means has opened up a fascinating domain of re- 
search [H- [i| Unlike in magnetic traps, where the spin 
of the atoms is effectively frozen, in an optical dipole 
trap the magnetic degrees of freedom dramatically af- 
fect the structure of the condensates. In particular, their 
equilibrium states exhibit richer degeneracy of physically 
distinguishable states than magnetically trapped BECs, 
with the degeneracy, or order, parameter depending on 
the atomic spin. In this paper we study a ferromag- 
netic spin-1 BEC in 87 Rb |fj whose degeneracyparame- 
ter is determined by the rotations of the spin (jf la] . We 
show that in the classical mean-field, or Gross-Pitaevskii, 
approximation there can exist a vortex line with a free 
end, terminating on a hedgehog-like spin configuration 
whose superfluid velocity profile coincides with the vec- 
tor potential of the Dirac magnetic monopole 0, [lnL Hll | . 
Magnetic monopoles are important in quantum cosmol- 
ogy; explaining their observed low density stimulated the 
theory of cosmological inflation [l^. The terminating 
vortex line, known as the Dirac string, is possible in a 
ferromagnetic spinor BEC because the noncommutativ- 
ity of the spatial rotations results in ambiguity of the 
condensate phase difference between two spatial points. 
We also show that a monopole can be attached to an 
anti-monopole with a vortex line, such that neither end 
reaches the boundary of the atomic cloud, forming a 
"dipole", Fig. HI 

In a magnetically trapped BEC the order parameter, 
determined by the complex scalar field, has a well-defined 
phase. Hence, its circulation around any closed loop can 



only equal an integer multiple of 27r, where the integer 
coefficient represents a topologically invariant winding 
number. As a result, any vortex line is either closed or 
terminates at the boundary of the atomic cloud. In a 
ferromagnetic spinor BEC the order parameter symme- 
try group is determined by the spatial rotations of the 
spin to be SO (3). This can be represented by a local co- 
ordinate axis, or a triad. A vortex line in a ferromagnetic 
BEC is determined by an integer multiple of 2tt spatial 
rotations of the triad about the direction of the spin, 
when one moves around any closed loop surrounding the 
vortex core. Due to the noncommutativity of the spatial 
rotations, such a vortex is not a topological invariant, 
but can be continuously deformed, e.g., to a vortex line 
with the opposite sign or to a disgyration. Since the rel- 
ative condensate phase between two spatial locations is 
not uniquely defined, a vortex line may also terminate in 
the middle of the atomic cloud. This is analogous to the 
possibility of vortex lines with free ends in the angular 
momentum texture of superfluid liquid 3 He-A 

In this paper we consider the mean-field model of a 
ferromagnetic spin-1 BEC. We numerically integrate the 
Gross-Pitaevskii equations in real time to find the dy- 
namics, and in imaginary time to find energetic minima. 
Our initial states are based on analytic expressions for the 



FIG. 1: Spin fields (F) for the monopole of Eq. @) (left) and 
the dipole of Eq. J7J (right). In all figures the positive z axis 
is upward, and cones are proportional to (F). The dipole's 
singular points at c = ±(2.5/V2)xh are near the second and 
fourth planes of spin vectors. 
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monopoles and dipoles, with overall Thomas-Fermi den- 
sity profiles. In 3 He studies dipoles have been referred 
to as monopolium 01 . In that context it was suggested 
that monopolium lattices might be dynamically stable in 
rotating systems. 

Previous theoretical work on monopoles in atomic 
BECs has been limited to anti- ferroma gne tic systems 
[lil Hil H?| . two-component monopoles |l8|. or analo- 
gous two-dimensional structures called "spin monopoles" 
|l9l |20| . This work has shown that monopoles may be 
created using phase im p rintin g techniques for creating 
solitons and vortices [Tj| W& 123 • Spin textures have 
also been recently experimentally created in an atomic 
spin-1 23 Na BEC by means of rotating the spin profile 
with spatially dependent external magnetic fields |23| . 



MEAN-FIELD MODEL FOR MONOPOLES AND 
DIPOLES 

The second quantized Hamiltonian for the spin-1 BEC 
in an optical dipole trap reads |7Hi. 153 . E l Eg 



H = I d 3 r - — • V^ Q + V(r)^J a 
[ Ira 

+ ^ c 2^a^ F a 7 ' F PSM 7 \ , (1) 



where ip a is the field annihilation operator for the a 
Zeeman component and m is the atomic mass. V(r) 
is the optical trapping potential. In experiments, op- 
tical traps with a wide range of aspect ratios and depths 
have been demonstrated |27j. In this paper we consider 
an isotropic harmonic trapping potential with frequency 
to: V(r) = muj 2 r 2 /2. The magnetic field vector is de- 
noted by B, [Ib is the Bohr magneton, and the Lande 
factor gp = —1/2 for the F = 1 transition of 87 Rb. In 
Eq. , F is the vector formed by the three components 
of the 3x3 Pauli spin-1 matrices [24|, and we have used 
the Einstein summation convention over repeated indices. 
The Zeeman energy is correct to first order in the mag- 
netic field. Co and C2 are the spin-independent and spin- 
dependent two-body interaction coefficients. In terms of 
the s-wave scattering lengths clq and 02 for the channels 
with total angular momentum zero and two they are: 
Co = 4:Trh 2 (ao + 2a2)/3m and C2 = Airh 2 (a,2 — ao)/3m. 
We use the values given by van Kempen et al. [2^] 
for 87 Rb: ao = 101. 805 and = 100. 4as, so that 
(a + 2a 2 )/3 = 100.9a B and (a 2 - a )/3 = -0.47a s , 
where the Bohr radius as — 0.0529 nm. We also use 
uj = 2tt x 10 s -1 , and assume a total number of N = 10 6 
atoms. We make our lengths dimensionless with the 
length scale x\ w = \J h/muj. As usual, the dynamics 



is invariant up to changes of length and time scale under 
changes of u> and N provided N 2 ui is constant. 

In the mean-field Gross-Pitaevskii approximation the 
order parameter is the three-component spinor \I/(r) = 
^/n(r)C(r), with C(r) t C( r ) = E and n(r) is the total 
density. Then the expectation value of the spin is given 
by (F) = QFapCp- The contribution of the two-body 
interaction to the energy is 
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d 3 r\-n 2 (c + c 2 (F) 
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The preceding parameters for 87 Rb correspond to the 
ferromagnetic case, since C2 < 0, and the energy is min- 
imized when I (F) | = 1 throughout the BEC for the case 
of a uniform spin distribution. 

For the ferromagnetic state |(F)| — 1, all the degener- 
ate, but physically distinguishable, states are related by 
spatial rotations of the atomic spin. These may be con- 
veniently defined by means of the three Euler angles : 

C = I Co I =e iv e- iF * a e- iF «l 3 e- iF * T | 
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where ip' = tp — r and ip denotes the macroscopic conden- 
sate phase. The combination p — r indicates an invari- 
ance under the change of the macroscopic phase if it is 
simultaneously accompanied by a spin rotation r of the 
same magnitude. The equivalence of the phase change 
and the spin rotation represents a broken relative gauge- 
spin symmetry. Consequently, the symmetry group of 
the order parameter is fully determined by the Euler an- 
gles to be SO (3). The direction of the expectation value 
of the spin is (F) = zcos/3+sin/3(xcosa-|-ysina), where 
x,y,z are the Cartesian coordinate unit vectors. 

We are interested in the ferromagnetic state with a sin- 
gular pointlike core in the atomic spin field. Using the ra- 
dial, azimuthal, and polar spherical coordinates (r, 0, </>), 
we select j3 = 0, a = <f», and <p' = —<fr m Eq. ©, yielding 



c = 



"^cos 2 (0/2)\ 
'sin (9)/V2 
in 2 (9/2) j 



(4) 



The corresponding spin exhibits a radially outward 
hedgehog field and forms the monopole (F) = z cos 9 + 
sin 9(x cos <fi+y sin 4>) . However, for 9 = the value of <j> is 
not defined and the first component is singular. There- 
fore, the point singularity, or the hedgehog, at r = 
is attached to a line singularity (in this case a doubly- 
quantized SO (3) vortex) along the positive z axis. This 
is called a Dirac string. Note that for z < there are 



3 




the dipole spinor 



FIG. 2: Spinor densities for the monopole (left) and dipole 
(right). The spinors are given by Eq. Q and Eq. J7J respec- 
tively, multiplied by the Thomas-Fermi and vortex density 
profiles, as described in the text. All lobes are bounded by 
density isosurfaces for a spinor component for x < 0. On the 
x = plane the isosurfaces are capped with a density col- 
ormap for the corresponding spinor component. Monopole: 
the top lobe is for |*l/i| 2 . The density zero associated with 
the doubly quantized vortex is seen along the positive z axis. 

The bottom lobe is for | \l/ i | 2 , and the central lobe is for 

l^o] 2 - The singly quantized vortex can be seen. Dipole: the 
top and bottom lobes are for l^i | 2 , the central one for l^-i | 2 , 
and the remaining two lobes for l^ol 2 - The vortices in vp-i 
and can be seen on the z-axis. 



no line singularities. The structure of the monopole in 
terms of the individual spin components is: a doubly- 
quantized vortex line in the \&i component, a vortex with 
unit circulation in the component, and no vortex in 
the ^-i component. The corresponding spinor densities 
are shown in Fig. 

The corresponding superfluid velocity is 



M 



M 



(Vy/ - cos/3Va) 



Mr sin 9 



(5) 



where </> is the polar unit vector. This has the form of 
the electromagnetic vector potential of the Dirac mag- 
netic monopole, analogous to the corresponding angular 
momentum texture in superfluid liquid helium-3 13]. 

An anti-monopole with spin vectors pointing radially 
inward (F) = — zcos# — sin#(icos</> + j/sin</>), and the 
Dirac string along the negative z axis, is given by 



-2i<j> gin 2 

— i<f> 



(0/2) \ 
sin(0)/\/2 
; 2 (0/2) J 



(0) 



An interesting system results from attaching a 
monopole to an antimonopole by joining their Dirac 
strings. The resulting vortex line with two free ends then 
no longer reaches the boundary of the atomic gas. From 
the spin field given by Soni |14| follows an expression for 




l + (r 2 cos2(9-c 2 )/A 
^V2(r 2 sin26>)/A 
- (r 2 cos2(9-c 2 )/A}, 



where 



A 



(J) 



(8) 



This has monopole and antimonopole type singularities 
on the z-axis at z = \c\ and z = — |c| respectively. On the 
z-axis for \z\ > \c\, (i = 1 and C-i = 0, while between 
the poles, i.e. for \z\ < |c|, £i = 0, and £-i is singular, 
representing a doubly-quantized SO (3) vortex line. The 
corresponding spinor densities are shown in Fig. 

The topological stability of linear and pointlike topo- 
logical defects may be characterized by the first and sec- 
ond homotopy groups. Since the second homotopy group 
of SO (3) is trivial, there are no stable topological point 
defects. However, as we previously showed, we may form 
the monopole by attaching it to a singular vortex line. It 
is not possible to form an isolated point defect in SO (3). 
The 'charge' W of the Dirac monopole is similar to the 
topological invariant of the point defects, which is de- 
termined by the way the order parameter behaves on a 
closed surface enclosing the defect 
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d(F) d(F) 



dxj dxk 



(9) 



Here e^-fe denotes the totally antisymmetric tensor. The 
integration is defined over a closed surface enclosing the 
monopole. The singular vortex line always cuts through 
the surface. However, as we show in the following section, 
the Dirac monopole can decay to a nonsingular coreless 
texture in which case (F) is well-defined everywhere. 



NUMERICAL RESULTS 

We now turn to the numerical studies of the clas- 
sical mean-field theory for the ferromagnetic spinor 
BEC. We integrate the relaxation and dynamics of the 
monopole and the monopole-antimonopole pair. The 
Gross-Pitaevskii e quat ions for the dynamics of the spinor 
mean-field are [2II2H 



*«* 7 . (10) 

The explicit forms of the two-body interaction terms, for 
a = 1, 0, —1 are, respectively, 

(c + c 2 )n<Fi - 2c 2 |*_i| 2 *i + c 2 <T ^ , 
(c + c 2 )n* - c 2 |* | 2 *o + 2c 2 *5* 1 *_ 1 , 



(c + c 2 )n*_i - 2c 2 |*i| 2 *_i + c 2 *^ 



■ 



(11) 
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We solve these equations numerically in three spatial di- 
mensions starting from initial monopoles and dipoles to 
find the dynamical evolution of the system. We also solve 
them in imaginary time t = —it to find minima, cither 
local or global, of the system's energy functional. This 
requires normalization of the Gross-Pitaevskii spinor af- 
ter each time step. Our numerical method is pseudo- 
spectral, using FFTs to evaluate the Laplacian, and a 
forth-order R unge -Kutta time step. It is described in 
detail in Ref. |29| . All our simulations were performed 
using 128 3 spatial grids for each of the three spinor com- 
ponents. The physical region simulated was centered on 
the trap, and was typically (25/v2)#/ia long in each di- 
mension. A typical time step was 10 _3 /c<j. 

The initial conditions for our numerical integrations 
of the Gross-Pitaevskii equations IjlUI) are either the 
monopole Eq. or dipole Eq. (0) with an overall 
Thomas- Fermi density profile n(r), multiplied by an ap- 
proximation to the vortex density profiles |24j, see Fig.0 
In both real and imaginary time the structures initially 
evolve by the filling in of the vortex core on the z axis by 
the vortex- free component: for the monopole this is the 
component, and for the dipole the Vfi component. 

The imaginary time evolution relaxes towards the low- 
est energy component. In the absence of a magnetic field 
this is the vortex-free component; ^-i for the monopole. 
This represents the relaxation of the monopole configu- 
ration towards a uniform spin texture. A magnetic field 
along the z axis can increase the energy of the vortex-free 
spin component, reversing the direction of the relaxed 
spin. However, the vortex cores may still be filled by 
the vortex-free component VP-i. This results in a core- 
less vortex, similar to the Anderson- Toulouse (AT) spin 
texture in superfluid liquid helium-3 0, 130L l3ll |. How- 
ever, the system no longer remained in the ferromagnetic 
ground state and regions of reduced spin magnitude, close 
to | (F) | = 0, were found in both the imaginary and real 
time simulations, see Figs. 01 and 0] 

The relaxation of the singular monopole to a nonsin- 
gular coreless texture may seem surprising at first, but 
it is topologically allowed. We can see how the singular 
monopole may be deformed to a nonsingular spin tex- 
ture by continuously rotating the direction of the spin at 
the positive z axis. Consider the following parametriza- 
tion: 

/e-^cos 2 [f3(p,z,s)/2}\ 
(= e-«Bm\p(p,z,s)]/y/2] , (12) 
V sin 2 [/3(p,z, S )/2] / 

where (p, <fi, z) are the cylindrical coordinates. We write 
for z > 0: (3 — arctan(p/z) + sir exp (— Cp/z), where 
C is here chosen to be much larger than the inverse of 
the radius of the atomic cloud. We continuously evolve 
the parameter < s < 1. At s — 0, we have the Dirac 
monopole with the singularity at p — 0, z > 0. The 
singularity can be removed by continuously evolving s 
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FIG. 3: Spin magnitude |(F)| and spin field (F) of the relaxed 
state, with a magnetic field B = \/21uo / p,B , of either an ini- 
tial monopole (B in z direction) or dipole (B in —& direction). 
The imaginary time simulations were run for —it — 30/u>, 
with a time step 0.002/ti;. Other parameters are as given in 
the text. Left. 3D slice plot of |(F}| in the central region. 
The width shown in each dimension is half of the total sim- 
ulation width of (25/y/2)xho- The slice planes are x = 0, 
y — 0, and z = 0. The colormap corresponds to |(F}| ranging 
between 0.046 and 1. Right. (F), for the dipole, on the plane 
z = 0. Quarter the simulation width is shown in the x and y 
dimensions. 

from to 1. At s = 1, the spin has rotated by the value 
of 7T at p — 0, z > 0, and the singularity has disappeared. 

Depending on the functional dependence of the param- 
eter (3 we may have different coreless textures. If (3(p) 
depends on p alone, with (3(0) = ir and (3 = at the 
boundary of the atomic cloud, we have the AT coreless 
vortex. On the other hand, with (3 = ir/2 at the bound- 
ary, we obtain the Mermin-Ho texture [32( . For the core- 
less textures, W, in Eq. ©, would represent a topolog- 
ically invariant winding number, provided that (F) was 
asymptotically restricted on a plane. However, unlike 
in superfluid liquid helium-3, where the interactions be- 
tween the superfluid and the walls of the container may 
fix the boundary values of the spin, in atomic BECs the 
spin may generally rotate and change the value of W. 

The continuous deformation of the singular Dirac 
monopole to a nonsingular texture is related to the topol- 
ogy of line defects in ferromagnetic spin-1 BECs. The 
ground state manifold, with the symmetry group SO (3), 
corresponds topologically to a three-sphere (a sphere in 
four dimensions) with the diametrically opposite points 
set to be equal, S 3 /Z 2 - In this space the only closed 
paths which cannot be continuously deformed to a point 
are precisely the paths joining the diametrically opposite 
points of the three-sphere |3jj. Hence, only vortices with 
winding number equal to one are topologically stable 
and, e.g., a vortex line with winding number equal to two 
can be continuously deformed to a vortex-free state. 

The regions of low spin in the AT vortex represent the 
largest values of the spin gradient in the texture, indi- 
cating a strong order parameter bending energy. The 
reduced spin value results from the very weak ferromag- 
netic energy of the 87 Rb spin-1 condensate: it is ener- 
getically more favorable for the system to violate the fer- 
romagnetic spin constraint than to create a spin texture 
with a very large spin gradient. The regions with a large 
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variation in the spin orientation continuously mix the fer- 
romagnetic |(F) | = 1 and the polar |(F)| = phases of 
the BEC, see Fig. [21 While the violation of the spinor 
phase and the variation of the spin value has been previ- 
ously investigated in the case of singular defects , it 
is interesting to observe that this may happen even as a 
result of the order parameter bending energy of nonsin- 
gular textures. 

In the presence of cylindrically asymmetric perturba- 
tions wc found the coreless AT texture to be unstable 
with respect to drifting towards the boundary of the 
atomic cloud. This is analogous to the slow drift of a 
vortex line in a single-component BEC in a nonrotating 
trap 0. We did not perform any quantitative studies 
to determine whether the drift of the spin texture could 
possibly be suppressed by the external magnetic fields. 

The dipole relaxes to a state which is the same as in 
the case of the monopole, apart from reversed spin di- 
rections due to the reversed magnetic field. Although 
the monopole and dipole relax to the same type of state, 
their imaginary time evolution follows different paths. In 
particular the cylindrical shell of low spin magnitude de- 
velops differently. For the monopole the transition from 
(F) = — z in the core to (F) = z at the boundary occurs 
by the spin vectors in the interfacial region rotating ra- 
dially outward. For the dipole, however, the direction of 
spin rotation differs in the regions z > and z < 0. The 
spin rotates radially outward for z > 0, and radially in- 
ward for z < 0. To accommodate these different rotation 
directions, near z = the spin magnitude drops close to 
zero early in the relaxation. The cylindrical shell of low 
spin magnitude then grows in the ±z directions as the 
relaxation progresses. In contrast, for the monopole the 
growth occurs from the gas boundaries towards the z — 
plane. 

The real time evolution of the dipole and monopole, 
without any applied magnetic field, are shown in Figs. 01 
and[3 The initial states were first relaxed for an imagi- 
nary time it = 0.1/w before the real time evolution. This 
ensured that the total density was close to the steady 
state distribution, hence suppressing density oscillations. 
The dynamics is quite complex and hence we restrict our- 
selves to describing the major features. The spin field of 
the dipole is generally more robust in the dynamical evo- 
lution than that of the monopole. Apart from the filling 
of the vortex core by the vortex- free component, recog- 
nizable structures survive around a half period (t = tt/lu) 
of evolution. However, the original states are noticeably 
different after about a period of evolution, which for our 
parameters is 2ir/u> = 0.1 s. This is likely to be an upper 
bound on the time available to experiment with these 
structures. At tot = 18, approximately three harmonic 
oscillator periods, the spin field of the initial monopole 
has developed an approximate dipolar structure, Fig. 
However the spinor is completely different to the dipole 
of Eq. Q. A major difference between the dynamics of 
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FIG. 4: Spin fields (F) (left column) and spin magnitudes 
|(F) | (right column) after real time evolution of the dipole of 
Fig- Q for uit = : 6 (top row), and 12 (bottom row). The 
spin magnitudes are shown on slice planes x = 0, y — 0, and 
z — 0. The colormap corresponds to |(F)| ranging between 
and 1. At ut = 6 the minimum of |(F)| is 0.0033. At ut = 12 
it is 0.0019. Approximately the central three-quarters of each 
simulation dimension is shown. Parameters given in text. 
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FIG. 5: Spin fields (F) after real time evolution of the 
monopole of Fig. Q for uot = : 3 (top left), 6 (top right), 12 
(bottom left), 18 (bottom right). Parameters given in text. 



the monopole and dipole is that the dipole develops large 
regions centered around the z = plane where the spin 
magnitude |(F)| drops to zero, see Fig. QJ 

We have examined the relaxation and dynamics of 
Dirac monopoles and dipoles in trapped ferromagnetic 
spinor BECs, for realistic parameters of 87 Rb. There 
are several interesting aspects to be explored further. 
We have not investigated the dynamical stability of the 
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monopole-antimonopole pair or the effects of rotation. As 
proposed by Soni [l4|, analogous structures, although not 
energetically stable in the presence of dissipation, could 
still be dynamically stable in rotating superfluid liquid 
hclium-3. In the absence of dissipation, our dynamical 
evolution also suggests the dipole spin field to be surpris- 
ingly robust. However, as our studies have shown, due 
to the weak ferromagnetic energy of 87 Rb atom, any dy- 
namically stable dipole configuration would also include 
regions of significantly reduced spin magnitude, present- 
ing a unique and highly nontrivial topological object for 
further studies. 
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